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Abstract 

The van der Waals (VDW) equation of state describes a thermal equilibrium in system of particles, 
where both repulsive and attractive interactions between them are included. This equation predicts an 
existence of the 1st order liquid-gas phase transition and the critical point. The standard form of the 
VDW equation is given by the pressure function in the canonical ensemble (CE) with a fixed number 
of particles. In the present paper the VDW equation is transformed to the grand canonical ensemble 
(GCE). We argue that this procedure can be useful for new physical applications. Particularly, the 
fluctuations of number of particles, which are absent in the CE, can be studied in the GCE. For the 
VDW equation of state in the GCE the particle number fluctuations are calculated for the whole 
phase diagram, both outside and inside the liquid-gas mixed phase region. It is shown that the scaled 
variance of these fluctuations remains finite within the mixed phase and goes to infinity at the critical 
point. The GCE formulation of the VDW equation of state can be also an important step for its 
application to a statistical description of hadronic systems, where numbers of different particle species 
are usually not conserved. 

PACS numbers: 12.40.-y, 12.40.Ee 
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I. INTRODUCTION 


Statistical approach appears to be rather successful in a description of hadron multiplicities 


in high energy collisions. It leads to surprisingly good agreement between the results of the 
hadron resonance gas and the experimental data (see, e.g., Refs. M) . The excluded volume 
(EV) model introduces the effects of hadron repulsion at short distances. The EV procedure 
substitutes a system volume V by the available volume 14v = V — bN , where b is the particle 
proper volume parameter and N is the number of particles. This substitution can be easily 
done in the canonical ensemble (CE) formulation, where the (conserved) particle number N is 
an independent variable. The ideal gas pressure Pid is then transformed to the EV pressure. 


Pid(R,T,iV) = 


NT 

'V 


^ PeAV,N,T) = 


NT 


V - bN ' 


( 1 ) 


w 


rich is the function of volume V, temperature T, and number of particles N. It was shown 

u n 

P, 18| that the substitution of V by 14v leads to a shift of the chemical potential p and to a 
transcendental equation for the pressure function in the grand canonical ensemble (GCE): 


Pev(T,/i) = Pid[T,/i -pev(R,T,/r)] , 


( 2 ) 


where Pid(T, p) is the pressure of the ideal gas calculated in the GCE. The EV model ([2]) can 
be equivalently formulated in terms of the mean-held (see Refs, js ill). This makes it possible 
to incorporate other hadron interactions within unihed mean-held approach. The EV hadron 

n n 

equation of state is used for htting the hadron multiplicities (see, e.g.. Refs. |12l. 11311) and for 


the hydrodynamic model calculations of nucleus-nucleus collisions (see, e.g.. Refs. jl4j-ll6|) 


In the present paper t 


terms (see, e.g.. Refs. 


17 


le equation of state with both repulsive (6 > 0) and attractive (a > 0) 

3 ), 


p{V,T,N) = 


NT 


V - bN 


N^ 


( 3 ) 


is considered. The equation of state (E]) was suggested by van der Waals (VDW) in 1873, 
and for his work he obtained the Nobel Prize in physics in 1910. In the present paper the 
VDW equation ([3]) will be transformed to the GCE. The GCE formulation will be then used 
to calculate the particle number huctuations in the VDW gas. 
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Normally, the GCE is an easiest one from the technical point of view. Thns, it is nsed more 
freqnently, and a typical problem in statistical mechanics is to transform the GCE resnlts to 
the CE ones. However, the pressnre fnnctions of the EV ([T]) and VDW ([3]) have their simplest 
and transparent forms jnst in the CE. A comparison of a simple analytical form ([T]) with a 
transcendental eqnation ([2]) elncidates this point. On the other hand, the CE pressnre does not 
give a complete thermodynamical description of the system. This is becanse the parameters 
V, T, and N are not the natural variables for the pressnre fnnction. 

Does one really need the GCE formnlation for the statistical description of hadrons? There 
are at least three reasons for a positive answer to this qnestion. 

First, the CE pressnre, e.g., Eq. ([T]) or Eq. ([3]), does not give any information concerning the 
particle mass m and degeneracy factor d. However, the average energy of the system depends 
on the valne of m, e.g., the thermal part of energy per particle eqnals to m + 3T/2 at m S> T, 
and to 3T at m <^T. 

Second, the nnmber of hadrons of a given type is nsnally not conserved. For example, the 
nnmber of pions can not be considered as an independent variable and is a fnnction of volnme 
and temperatnre. The GCE formnlation of the VDW eqnation can be therefore an important 
step for inserting the repnlsive and attractive VDW interactions into the mnlti-component 
hadron gas. 

Third, Eqs. ([T]) and ([3]) of the CE give no answer abont the particle nnmber flnctnations. 
Formally, the nnmber of particles does not flnctnate in the CE: the total nnmber of particles 
Nq is constant in the fnll volnme Vq. However, the nnmber of particles N starts to flnctnate if 
one considers a snb-system with V < Vq. If V <C Vq the V-flnctnations shonld follow the GCE 
resnlts. For the EV model the particle nnmber flnctnations were calcnlated in Ref. . 

The paper is organized as follows. In Sec. [TTl we discnss how to transform the CE pressnre 
to the GCE. This procednre is illnstrated by a derivation of the EV model for pev(E,/i) in the 
GCE. In Sec. [TTIlthe VDW eqnation of state ([3]) is transformed to the transcendental eqnation 
for particle nnmber density n(T, /i) in the GCE. The particle nnmber flnctnations and their 
behavior in a vicinity of the critical point is then analyzed in Sec. IIVI A snmmary in Sec. |V] 
closes the article. 


3 


II. FROM CANONICAL ENSEMBLE TO THE GRAND CANONICAL ENSEMBLE 


In the GCE the macroscopic states of a thermal system are dehned in terms of the grand 
potential O which is a fnnction of 3 physical variables - volume V, temperature T, and chemical 
potential p. In a spatially homogeneous system the grand potential O is straightforwardly 


connected to the system pressure 


17| 


0(1/, T,p) = -p{T,p)V . 


(4) 


Therefore, pressure p as a function of its natural variables^ temperature T and chemical potential 
/r, gives a complete description of the equation of state. One calculates the entropy density s, 
(conserved) number density n, and energy density £ as 


4T,,) ^ 1 I 


n(T,/r) = 


dp\ 
dp,) T 


e{T,p) = Ts + pn — p . (5) 


The thermodynamical potential in the CE is the free energy F{V, T, N) (or Helmholtz poten¬ 
tial). The function F depends on its natural variables, volume V, temperature T, and number 
of particles N, and it gives a complete thermodynamical description of the system properties. 


^jU 

m 


ropy S, pressure p, chemical potential p, and total energy E are calculated in the CE as 


s(r,T.]V) = -(^l 


V,N 


p(r,T.]V) = -(|/ 


T,N 




and 


E{y,T,N) = F + TS . (7) 

If one knows p = p{V, T, N) in the CE it is not yet possible to hnd S{V, T, N), p{V, T, N), 
and E{y,T,N) in unique way. For a complete thermodynamical description in the CE one 
needs to reconstruct free energy F, thus, the following equation should be solved: 

The solution of Eq. ([8]) can be presented as 

F{V,T,N) = F{Vo,T,N) - f dV'p{V',T,N) . (9) 

dVo 


4 





To find F{V,T, N) one needs to know not only p{V,T,N) but, additionally, the free energy 
F{Vq,T, N) as the function of T and N at some hxed value of Id = Vq. We assume that for 
hxed T the role of particle interactions becomes negligible at very small density N/Vq 0. 
Thus, the ideal gas expression, F{Vo,T, N) = T , N), becomes valid at Vq ^ cxd. For 


relativistic ideal Boltzmann gas the free energy reads 


Fid(l/,r,7V) = - NT 


1 + In 


li 


N 


where 


(l){T;d,m) = 


d , 2 ,, / / 7 T 5 - o driFT /m 

— kdk exp(-VP + mVT) = (- 


( 10 ) 


( 11 ) 


In Eqs. flTU]) and flTT]) . d is the degeneracy factor, m is the particle mass, and quantity (j) has a 
physical meaning of the GCE particle number density at zero chemical potential. 

Let us start with an example of the EV volume equation ([1]) in the CE. The particle proper 
volume parameter in Eq. ([1]) equals to 6 = 4 • (dvrr^/S) with r being the corresponding hard 
sphere radius of particle. This equation describes the repulsion between particles and can be 
rigorously obtained (in particular, a factor of 4 in the expression for b) for a gas of the hard 
balls at low density (see, e.g.. Ref. jf^). 


From Eq. 


it follows 


F{V, T,N) = - NT 


1 + In 


Vo(j){T; d,m) 


= - NT 


1 + In 


Vo4>{T;d,m) 

N 


N 


+ NT In 


Rn 


-0 

Jv V' - bN 

bN 


= Fid(R-6iV,T,iV) , 


( 12 ) 


V -bN 

where we have used that (Vq —6iV)/Vo —)■ 1 at Vq ^ cxo. For the free energy flT^ the chemical 
potential can be calculated as 


= 


fdF\ 


\dN ) 


V,T 


(V-bN)4,{T-,d,m) , ^ NT 

= -T 1 “-^. 


(13) 


For particle number density n = N/V one hnds from Eq. fll3p for the ideal gas, i.e., at 6 = 0: 

y = nid(r,/x) = exp(^^) (p{T;d,m) . (14) 

At 6 > 0 one can easily hnd from Eq. flT^ : 


^ — 'n^eviT, /i) 


nid(r, ii* 


1 + 6nid(T,/i*) 


fi = fi 


Ue^jT, /i)T 

- 6nev(T,/i) 


(15) 
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or equivalently, 


n (T u)T 

Pev{T,fJ.) = Pid(T,/i*) , /i* = /i - b- - - bpeviT,fx) . (16) 

1 - bne^{T,p) 

Therefore, to solve the excluded-volume problem in case of the Boltzmann statistics one has 
to obtain solution of transcendental equation ffTHjl with respect to particle density nev{T,p) or 
of transcendental equation ffT6|) with respect to pressure PeviT,fi) for every point in (p, T)- 
plane, i.e^ in the GCE. These equations coincide with those suggested in Refs. BJni and in 


Refs. 


71, Isj, respectively. 


III. VAN DER WAALS EQUATION OF STATE IN THE GCE 


GCE formulation 


The VDW equation for the system pressure is dehned in the CE in the form of Eq. 
Using Eq. ([9]) one obtains 


fVo 

F(U,T,V) = Fid(Uo,r,7V) + / dV' 

Jv 


NT 


V'-bN 


'v^ 




^ Fid(U-&iV,T,iV) - a— . 


Chemical potential then reads 
h = 

and from Eq. flT8|l one obtains 


dF \ 

mJvT 




N 


N 


V-bN U 


N _nid(T,p*) 

T/ - 1 p) -I I 7 /rp 7 

V 1 + 6nid(T,p*) 


p* = p — b 


n T 
1 — bn 


+ 2a n . 


(17) 


(18) 


(19) 


The function n(T, p) which is the solution of transcendental equation flT^ is reduced to ne^(T, p) 
flTSl) . if a = 0, and to ni^iT^p) ffTT|) . if a = 6 = 0. At hxed T and p, a presence of parameter 
b > 0, which describes the repulsion between particles, leads to a suppression of particle number 
density n{T,p), while attractive interactions described by a > 0 lead to its enhancement. 
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The energy density of the VDW gas can be calculated within the GCE from the pressure 
function, p{T, /i), 


£(r./<) = r 


dp 

+ - p = 


£id(r,p' 


1 + 6nid(T,/i*) 


— an 


an 


n 


= Cid (T;m) 

where J1 is defined in Eq. flT^ and 

+ /c^ exp \/w? + k'^/T) 

6id(r,T/l) — 


= 3T + m 


Ki{m/T) 


( 20 ) 


( 21 ) 


k'^dk exp (— + k'^/T) K 2 {ni/T) 

is the average energy per particle in the ideal gas. A comparison of the right hand side of 
Eq. fl2U]) with the corresponding ideal gas expression, = eid(T; m) nid, demonstrates the role 
of the a and b parameters for the system energy density. Both a and b parameters influence the 
particle number density and transform njd into n: a > 0 leads to an enhancement of particle 
number density and 6 > 0 to its suppression. Parameter b does not contribute however to the 
average energy per particle, whereas the mean field —an created due to parameter a is added 
to the average energy per particle and changes it from Cid to e = Cid — a n. 

It is also instructive to consider another derivation of Eq. fll9p which based on the procedure 


used in Ref. 


20|. Introducing particle number density, n = N/V, the VDW pressure function 


can be written as 


p{T,n) = 


n T 


— an 


( 22 ) 


1 — bn 

For the particle number density in the GCE the following differential equation can be obtained 


n{T,p) = 


( dp{T,p) 

V dp 

T 


[ dp{T,n)\ [ dniT^p) 


T 

— 2a n 


\ dn 


(1 — bnY 

A general solution of Eq. fl2^ can be rewritten as 

T 


( dn{T,l) 

V dp 




dp 


p = dn 


n{l — bnY 


2a 


= T In 


n 


+ 


bT n 


2a n . 


(23) 


(24) 


(j){T;d,m) {1 — bn) \ ' 1 — bn 
The temperature dependent integration constant in the solution fl2T]) is dehned in the form 
which guarantees the ideal gas solution flTT)) at 6 = a = 0. One can easily check that functions 
n{T,p) dehned by Eq. fIMl) and Eq. flT^ are identical. 
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Figure 1: (Color online) (a): The dimensionless form of the VDW isotherms for pressure, p, versus 
the volume per particle v = n~^. The dashed-dotted lines present the metastable parts of the VDW 
isotherms at T < 1, whereas the dotted lines correspond to unstable parts. The full circle on the 
T = 1 isotherm corresponds to the critical point, (b): The phase diagram for the VDW equation of 
state on the (n, T) plane. The phase coexistence region resulting from the Maxwell construction is 
depicted by grey shaded area. 


B. Phase Diagram for the VdW Gas 


The VDW equation of state contains a 1®* order liquid-gas phase transition and has a critical 
point. The thermodynamical quantities at the critical point are equal to Q: 


T, = 


8a 


nr_ = 


3b 


Pc = 


2762 


(25) 


The VDW equation ([3]) can be then rewritten in the following dimensionless (reduced) form: 


8T n 


P = 


3n2 


n 


(26) 


where n = n/nc, p = p/pc, and T = T/Tc. In the dimensionless presentation fl26|) the VDW 
equation has a universal form independent of the values of a and b, and the critical point fl25|l 
is transformed to Tc = Pc = ric = I ■ 


The dimensionless VDW isotherms are presented in Fig. [T] (a) as functions oiv = n~^ . To 
describe the phase coexistence region below the critical temperature the VDW isotherms should 
be corrected by the well-known Maxwell construction of equal areas. These corrected parts of 
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the VDW isotherms are shown by the solid horizontal lines. Figure [T](b) depicts the liquid-gas 
coexistence region on the (n,T) plane. The pure gaseous phase exists at all T < 1. However, it 
can be hardly seen at T -C 1 in Fig. [1] (b). This is because of very small admittable values of n 
inside the pure gaseous phase at T —)■ 0, e.g., the largest value of the particle number density 
in the gaseous phase at T = 0.25 is n = 5 • 10“®. At any value of T < 1, the particle number 
density in the mixed phase is given by 

n = ^Ug + {1- ^)ni , (27) 

where ^ and 1 — ^ are the volume fractions of the gaseous and liquid components, respectively. 
The values of Ug and ni in Eq. fl27)l are the particle densities of, respectively, the gaseous 
and liquid phases at the corresponding boundaries with the mixed phase. For example, at 
T = 0.25 one finds rig = 5 ■ 10“^ and ni = 2.75. For ^ = 0.95 Eq. fl2711 then gives n = 
0.95 • 5 ■ 10“® -I- 0.05 ■ 2.75 = 0.14. Therefore, at n = 0.14 and T = 0.25 the gaseous component 
occupies 95% of the total volume of the mixed state but gives a negligible contribution to the 
total number of particles. Note also that n = 3 is an upper bound for the dimensionless particle 
density. 

Let us consider the Maxwell rule in more details. At each T < 1 the VDW loops appear in 
the isotherm p(n). They are shown by the dashed lines in Fig. [1] (a). Each loop consists from 
metastable part, where dp/dv < 0, and unstable part, where dp/dv > 0. They are shown in 
Fig. [T] (a) by the dashed-dotted and dotted lines, respectively. According to the Maxwell rule, 
the loop between values of vi and Vg is substituted by the constant value po (solid horizontal 

line in Fig. [U (a)), and the following condition should be fulhlled: 

/•A _ _ _ _ _ 

dvp = Po {vg - Vi) . (28) 

Jvi 

It can be proven that the GCE formulation forbids the Maxwell loops and it leads auto¬ 
matically to the Maxwell construction. Let us assume that the system can be divided into two 
parts - gas and liquid - with 

N N, 

Ng + Ni = N , Vg + Vi = V , and ^ ^ . (29) 

Vg Vl 

The free energy of the system can be then presented as 


F{V,T,N) = F{Vg,T,Ng) + F(V,T,W), 


(30) 
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where the surface free energy is neglected in the thermodynamic limit K —)■ cxo. An equilibrium 
state corresponds to a minimum of the free energy. This leads to 

dF\ f dF 


dV„ 


= 0 , 


dN„ 


= 0 


(31) 


' 9 / T,Ng \^-^'9/T,Vg 

Relations in Eq. (ITT]) are equivalent, respectively, to 

p{T,Vg,Ng) = p{T,Vi,Ni) , KT,Vg,Ng) = p{T,Vi,Ni) . (32) 

The hrst condition in Eq. (13^ means that the Maxwell loop is transformed to the mixed phase 
region with equal pressures. And the second condition in Eq. fl52]) leads in a unique way to the 
Maxwell rule of equal areas. The equilibrium conditions of phase coexistence are then fulhlled: 
both phases have the same temperature, pressure and chemical potential in the coexistence 
region. These conditions are known as the Gibbs rules for a 1®* order phase transition 181. 


IV. PARTICLE NUMBER FLUCTUATIONS IN THE VDW GAS 


A. Particle Number Fluctuations in the Pure Phases 


The variance of the total particle number in the GCE can be calculated as 


Var[N] = (N^) - {Ny = T 


(a{N)\ 


= ri/(A 

\dfly rp 


(33) 


\ / T,v 

where symbol (...) denotes the GGE averaging, and n{T,fi) is the particle number density in 
the GGE. The scaled variance for the particle number fluctuations is then 


VarlN] T fdn\ 


It will be convenient to rewrite Eq. flT^ as 

n = (1 — hn) nid(T, p) exp 


hn 


exp 


2an 


Taking the derivative of Eq. fl35|l with respect to p one obtains 

n -1 


dn n 
^ “ T 


(1 — bny 


2 an 


(34) 


(35) 


(36) 
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Figure 2: (Color online) The lines of constant values of the scaled variance uj[N] are shown on the 
(n, T) phase diagram, outside of the mixed phase region. 


Therefore, the scaled variance fl5T|l equals to 


uj[N] = 


2 an 


-1 


(37) 


It is clearly seen from Eq. fl57)) that in the VDW gas the repulsive interactions suppress the 
particle number fluctuations, whereas the attractive interactions lead to their enhancement. 
Note that for a = 0 the value of Ci;[iV] fl37|l coincide with the result for the EV model obtained 


in Ref. 


19|. 


The scaled variance fl37|l expressed in terms n and T equals to 

-1 




9 [(3 - n)= 


n 

AT 


(38) 


In Fig. |2]the lines of constant values of a;[A^] are shown on the (n, T) phase diagram outside of 
the mixed phase region. At any hxed value of T, the particle number fluctuations fl38|l approach 
to those of the ideal gas, i.e. a;[A^] = 1, at n —)■ 0, and become small, u[N] 1, at n ^ 3. 
As it should be, the scaled variance fl55]) is always positive for all possible values of n and T 
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outside of the mixed phase region. At T —)■ 0 this is insured by simultaneous limiting behavior 
of n —>■ 0 in the pure gaseous phase and n —>■ 3 in the pure liquid phase. 

At the critical point (T = n = 1) the scaled variance of the particle number fluctuations 
diverges in the GCE. To study the behavior of a;[A^] in a vicinity of the critical point we 
introduce the quantities r = T — 1 -C 1 and p = n — 1 <^1. Expanding fl55]) at small r and p 
and keeping only the lowest orders of them one hnds 


n;[Ar] ^ 


4 

9 


r + -p^ + rp 
4 


(39) 


In particular, 


4 16 

u[N] = -T~^ at p = 0, and a;[A^] = Tpzat r = 0. (40) 

9 27 

Note that thermodynamical parameters T and n correspond to points outside the mixed phase 
region. Thus, in Eq. (I40|) . parameter r is positive, while p can be both positive and negative. 

As was mentioned above, the VDW equation of state permits the existence of metastable 
phases of super-heated liquid and super-cooled gas. These states are depicted by the dash- 
dotted lines on the VDW isotherms in Fig. [T] (a). In metastable phases the system is assumed 
to be uniform and, therefore, one can use Eq. fl38|l to calculate particle number fluctuations 
in these phases. In Fig. [3] the lines of constant values of Ci;[V] are shown on the (n,T) phase 
diagram for both stable and metastable pure phases, while the unstable region is depicted by the 
gray area. It is seen that the scaled variance remains hnite, and diverges only at the boundary 
between the metastable and unstable regions. We recall that at this boundary dp/dn = 0, 
where p is the dimensionless CE pressure fl26|) . One can easily show using Eqs. fl26|l and fl38|) 
that Ci;[V] —)■ cxo when dp/dn = 0. Note that metastable regions of the equation of state can 
be reached within fast non-equilibrium processes, whereas the unstable region is physically 
forbidden. Note that the thermodynamical relations are not fulhlled in the unstable region, 
e.g., nonphysical behavior with u[N] < 0 is found in this region. 


B. Particle Number Fluctuations in the Mixed Phase 

Let us consider a point (n, T) inside the mixed phase. The volume V is then divided into 
the two parts, Vg = occupied by the gaseous phase with particle number density Ug and 
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Figure 3: (Color online) The lines of constant values of the scaled variance a;[A^] are shown on the 
(n, T) phase diagram for both stable and metastable pure phases. The boundary between stable and 
metastable phases is depicted by the thick black line, and the unstable region is depicted by the gray 
area. 

Vi = (1 — occupied by the liquid phase with particle number density u/. The total number 
of particles N equals to Ng + Ni. The average value and scaled variance for the iV-distribution 
are 


(JV) = {N,) + (JV,) = C Kn, + (1 - On,] , (41) 

- W) - W )— — W )— - W )—' ' ’ 

In calculating the variations of Ng and Ni distributions inside the the mixed phase, one should 
take into account the fluctuations of the corresponding volumes Vg = and V) = (1 — ^)ld. 
In fact, this is the fluctuations of the parameter In the thermodynamic limit V —>■ oo, one 
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finds: 


(7V2) - (iV,)2 (AT,) 


{N) 

W) - (Ni)^ 


{N) 

(N,) 


{N) {N) 

(NgNi) - (TV,)(TV,) 




= - 2 


UgUl 


V [(e^) - (01 , 


(43) 

(44) 

(45) 


(N) n 

where Ci;^[TV,] and Ci;^[TV,] in Eqs. fl4^ and fl44)) correspond to the hxed value of and can be 
calculated using Eqs. fl33|l and fl34)l . One hnds 


a; [TV] = 


+ 


n [(1 — bUgY 
{Ug -riifv 


2ang 

~Y~ 


-1 


+ 


(1 - io)ni 


n 


(1 - hnif 


2ani 

~T~ 


-1 


n 


[«t - <01. 


where the equilibrium value {^) = is dehned as 

rii — n 


0 = 


ni 


Un 


(46) 


(47) 


To calculate the variance (^^) — (^)^ in Eq. fl46|l we can return to the CE formulation and 
use Eq. flHOl) rewriting it as the following 


F{V,T,N;0 = F{^V,T,Ng) + F [{1 - 0V,T, N^] . 


(48) 


Note that in the limit E —)■ cxo assumed in our study the particle number densities n, and n, 
in the GCE are identical to their CE values as a consequence of thermodynamical equivalence. 
The probability distribution 1E(0 is proportional to exp[—F(V, T, TV; ,^)/T]. The function F is 
given by Eq. (145]) and s a function of ^ it can be presented as a power series expansion over ^ — 
in a vicinity of the equilibrium value ^o- This gives the ^ normalized probability distribution 
[C is a normalization factor) in the form 


W{i) = C exp 




= Cexp 


{i - 

2(t2 


(49) 


where 


a 


= T 


fd^F\ 




-1 


?=Co 


T 


UgT 


T 

2an^ 


dPg 

dV„ 


+ 


+ 


niT 


dpi 

dVi 


-1 


2an‘f 


V Leo(l-K)' ^0 {l-io){l-bniy 1-^0, 


-1 


( 50 ) 
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Figure 4: (Color online) The lines of constant values of a;[A^] are shown on the (n,T) phase diagram 
inside the mixed phase region. 

Using the ^-distribution one finds 

r) - (0^ = 

and, finally, 

(1 -^o)ni 


(51) 


r 7> T] ^O'^g 

= 


2an„ 


1 -1 


+ 


n [(1 — bngY T 

{rig-nif 


+ 


2an] 


n -1 


n L(1 “ bniY T 


n 


2ani 


+ 
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Introducing dimensionless variables fl26|l one can rewrite Eq. fl52|l as 
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(53) 


In Fig. mthe lines of constant values of a;[iV] (1551) are shown on the {n,T) phase diagram 
inside the mixed phase. At small T, the particle number fluctuations inside the mixed phase 


15 














































0.4 


0.3 


^ 1 — 0.2 

0.1 

0.0 

0 2x10'^ 4x10'" 6x10'" 8x10'^ 1x10'" 

n 

Figure 5: (Color online) The lines of constant values of are shown at small temperatures and 
densities on the (n, T) phase diagram both inside and outside the mixed phase region. 

are small, u[N] -C 1, almost everywhere. This in because the main contribution to fl5^ comes 
from the liquid component of the mixed state, and the corresponding fluctuations inside the 
liquid phase are always small at small T. However, even at T —)■ 0, the scaled variance 
remains close to the ideal gas value, a;[iV] = 1, in a small region near the boundary with the 
gaseous phase. This region is not seen in Fig. H] because of very small n-values. A part of this 
region with 0.2 < T < 0.4 (both inside and outside of the mixed phase) is shown in Fig. [5l 
To study the behavior of fl53|) in a vicinity of the critical point we introduce 0 < Pg = 
1 — Ug <C 1, 0 < p; = fp — 1 1, and 0<t = l — T<^1. Att—)-0 one hnds Pg = pi = 2\/t 

and 

1 fi 

a; [AT] - - . (54) 

Finally, the scaled variance a;[A^] is shown in Fig. Elfor all possible values of n and T, both 
inside and outside of the mixed phase. It is interesting to compare Eq. for the scaled 
variance calculated inside the mixed phase and Eq. flMj) outside the mixed phase. These two 
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Figure 6: (Color online) The lines of constant values of uj[N] are shown on the phase diagram for 
0 < T < 3 and for all possible n values, both inside and outside the mixed phase region. 

equations give the same results at the phase boundary. At the boundary between the gaseous 
and mixed phase (i.e., .^o = 1) both the second and third terms in the right hand side of 
Eq. (I53|) become equal to zero, and the hrst term in (l5^ coincides with Eq. fl38D . Similarly, 
at the boundary between the liquid and mixed phase (i.e., .^o = 0) both the hrst and third 
terms in the right hand side of Eq. 135]) become equal to zero, and the second term in fl53p 
coincides with Eq. fl38|l . It might seem from Fig. [6] that at small temperatures huctuations 
have a discontinuity at the boundary between the liquid and mixed phase. It is, however, not 
the case: as mentioned above, when approaching the boundary from inside the mixed phase 
(■Co 0), the third term in the r.h.s. of Eq. flS5]l vanishes, but at small temperatures it vanishes 
rapidly and only in the very vicinity of the phase boundary. Consequently, this behavior is seen 
poorly on the scale used in Fig. |Hi 
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V. SUMMARY 


The VDW pressure function ([3]) has the simple and familiar form in the CE. However, there 
are some reasons for its reformulation in the GCE. First of all, the CE pressure does not give a 
complete thermodynamical description of the system. Besides, the GCE has clear advantages if 
one will try to extend the VDW model to multi-component hadron gas. The GCE formulation 
is also necessary to hnd the particle number fluctuations in the VDW gas. 

In the present paper we consider the procedure of a transformation of the CE pressure 
function to the GCE equation of state. This procedure has a form of the hrst order differential 
equation for the free energy F as a function of V. A general solution of this equation includes 
the boundary condition at V = Vq which is taken in the form of the ideal gas free energy in the 
limit of N/Vq —)■ 0. This adds to the GCE formulation an information of the particle masses 
and degeneracy factors which was absent in the CE pressure function. Following this procedure 
we transform the CE VDW pressure ([3]) to the transcendental equation flT^ for the particle 
number density in the GCE. Other thermodynamical functions of the VDW gas in the GCE 
can be then obtained in a unique way. 

The particle number fluctuations are calculated for the VDW gas in the GCE formulation. 
An analytical expression for the scaled variance of particle number fluctuations is derived in 
terms of dimensionless particle number density n = n/ric and temperature T = T jT^. It clearly 
demonstrates that in the VDW gas the repulsive interactions suppress the particle number 
fluctuations, whereas the attractive interactions lead to their enhancement. In the pure phases 
far away from the critical point the particle number fluctuations approach to those of the ideal 
gas, i.e., u[N] = 1, at n <C 1, and becomes small, a;[iV] -C 1, at n —)■ 3. The scaled variance 
a;[iV] increases in a vicinity of the critical point and diverges at the the critical point n = T = 1. 
Admitting a presence of metastable states, one observes that the scaled variance a;[iV] diverges 
also at the boundary between the metastable and unstable regions, where dp/dn = 0. 

The fluctuations of partial volumes Vg and Vi occupied by the gaseous and liquid components 
have influence on the particle number fluctuations in the mixed phase. However, the particle 
number fluctuations remain hnite in the mixed phase. They start to increase in a vicinity of 
the critical point and diverge at the critical point reached from inside the mixed phase. 
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We hope that the present study of the VDW equation of state within the GCE formulation 
can be useful for several applications in hadron physics. In particular, it looks promising to 
consider the nuclear matter with the VDW equation of state corrected by the Fermi statistics 
effects. It is also interesting to study the particle number fluctuations for hadronic systems in 
a vicinity of the critical point. 
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